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EVOLUTION  EQUATIONS  OF  C3I: 

CANNONICAL  FORMS  AND  THEIR  PROPERTIES 

Paul  H.  Moose 

In  this  paper,  five  different  two-species  non¬ 
linear  evolution  equations  are  described  that 
represent,  (a)  mixed  attrition  Lanchester  combat 
with  resupply  and  (b)  four  different  models  of 
Information  war.  Their  dynamical  properties  are 
analyzed  and  It  Is  shown  that  four  of  the  five  are 
environmentally  unstable.  The  meaning  of  this  Is 
Interpreted  for  03*!,  counter-O31,  and  Intel¬ 
ligence  operations.  For  the  environmentally  un¬ 
stable  systems  It  Is  shown  that  the  relationships 
between  the  parameters  In  the  model  are  critical  in 
establishing  their  mode  of  behavior,  whereas  in¬ 
itial  conditions  are  critical  In  determining  the 
actual  trajectories  of  evolution. 

In  several  earlier  papers,  we  have  suggested  the  use  of  non-linear  evolu¬ 
tion  equations  to  model  the  dynamics  of  the  average  Information  of  each  of  two 
opposing  sides  engaged  In  conflict.  Moose,  [1980],  developed  the  properties  of 
one  particular  version  of  a  model  for  counter-C3;  counter-C3  refers  to  "de¬ 
liberate  acts  to  distort  or  deny  Information  or  to  confuse  or  deceive  one's  op¬ 
position."  In  that  paper  several  alternative  counter-C3  models,  as  well  as  a 
model  for  Intelligence  activities,  were  suggested. 
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In  a  subsequent  paper.  Moose,  [1983],  developed  a  four-species  model  to 
account  for  Interactions  between  the  combat  forces  and  the  Information  on  each 
side.  Some  of  the  most  general  properties  of  that  model  exposed  the  fact  that, 
unlike  the  earlier  counter-C^  model,  a  four-species  model  Is  environmentally 
unstable.  That  Is,  for  some  ranges  of  the  parameters  of  the  model,  the  system 
Is  unstable,  whereas  for  others  It  Is  stable.  By  contrast,  an  environmentally 
stable  model  Is  stable  for  all  values  of  Its  parameters.  (Passive  electrical 
and  mechanical  networks  are  familiar  forms  of  environmentally  stable  systems, 
whereas  automatic  control  systems  are  environmentally  unstable  systems.) 

Evolution  equations  have  been  used  to  model  phenomena  in  the  social 
sciences  before.  Lanchester,  [1916],  Initiated  the  use  of  coupled  first-order 
differential  equations  to  model  the  attrition  of  military  forces  engaged  In 
combat.  Rlcharson,  [1960]  attempted  to  fit  the  growth  of  military  arms  In 
"arms  races"  to  linear  evolution  models.  In  a  recent  survey  article, 
Intrll Igator,  [1982],  cited  eleven  examples  of  the  use  of  differential  equa¬ 
tions  and  systems  of  differential  equations  to  account  for  various  political 
phenomena  related  to  conflict  resolution. 

Non-linear  evolution  equations  are  particularly  Interesting  because  they 
have  more  than  one  stationary  point.  This  enables  us  to  duplicate  mathemati¬ 
cally  the  unusual  behavior  of  many  complex  phenomena  that  seem  to  be  attracted 
for  a  period  of  time  toward  a  particular  condition,  then  In  an  apparently  arbi¬ 
trary  and  unpredictable  way  change  to  another  characteristic  condition.  This 
branch  of  mathematics  Is  now  referred  to  as  chaos  theory  and  the  characteristic 
types  of  system  behavior  are  called  "strange  attractors." 

In  this  paper,  we  shall  study  the  properties  of  five  different  two-species 
systems.  One  of  these  is  a  Lanchester  type  combat  force  attrition  model  with  a 
mixture  of  area  and  directed  fire;  the  relative  mix  will  be  made  dependent  on 
parameters  which  we  presume  are  related  to  the  "goodness"  of  the  C3j  system 
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of  each  side.  In  addition  we  shall  study  three  counter-C^  models.  Including 
the  model  of  Moose,  [1980].  These  models  are  Intended  to  capture  the  essence 
of  "Information  war."  The  fifth  model  was  suggested  by  Moose  [1980]  to  de¬ 
scribe  the  struggle  between  two  sides  for  "Intelligence  Information."  These 
five  models  are  shown  to  be  of  three  cannonlcal  types.  The  stationary  points 
and  stability  of  each  type  are  derived.  Types  "one"  and  "three"  are  environ¬ 
mentally  unstable  and,  as  shown  by  example,  are  very  rich  In  the  nature  of 
their  behavior.  Type  "two"  Is  environmentally  stable. 
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The  Equations  and  Their  Interpretation 


Coupled  two-species  evolution  equations  are  of  the  general  form 
X  -  -Fx(X,  Y,  v)  *  Vx  ) 

.  ~  (1) 

Y  *  -FJX,  Y,  w)  +  V 

j  ~  j  J 

Fx  and  Fy  are  attrition  functions.  They  depend  on  the  state  variables 
X  and  Y  and  on  a  set  of  parameters  u.  Vx  and  Vy  are  replenishment  terms. 
They,  as  well  as  the  u  may  be  constants  or  time  dependent.  If  the  parameters 
and  replenishments  are  constants,  the  system  Is  said  to  be  autonomous.  We 
shall  consider  only  autonomous  systems. 

Mixed  attrition  of  combat  forces  Is  characterized  by; 

1)  A  rate  of  loss  due  to  directed  fire  which  Is  proportional  to  the  size 
of  the  opposition's  force. 

2)  A  rate  of  loss  due  to  area  fire  which  depends  on  the  size  of  the  op¬ 
position's  force  and  of  the  density  of  one's  own  force. 

3)  A  natural  rate  of  loss  due  to  accident.  Illness,  dlsertlon,  etc. 

4)  Reinforcement  and  resupply  of  the  combat  force  In  the  field. 

In  order  to  account  for  these  conditions,  let 

X  »  -bxX  -  ay  Hy  Y  -  a  (1  -  u  )XY  +  Vx 

.  (2) 

Y  -  -byY  -  ax  yx  X  -  ax(l  -  ux)XY  +  Vy 

be  the  specific  characterization  of  (1).  The  parameters  bx  and  by  are  the 
natural  loss  rates,  oy  and  ax  are  related  to  the  rates  of  fire  and  prob¬ 
ability  of  kill  of  the  weapons  systems  In  use,  and  and  ux  describe  the 
fraction,  of  Vs  and  X/s  forces  respectively,  that  are  aimed  or  directed  fire. 
If  a  u  Is  zero,  then  all  the  fire  Is  area  fire.  If  a  u  Is  one,  then  all  the 


fire  Is  aimed.  The  Vx  and  Vy  are  the  average  rates  of  reinforcement  and 
resuppl  y  of  and  X« 

We  shall  be  Interested  In  analyzing  the  behavior  of  (2)  In  the  two-dimen¬ 
sional  parameter  space  wx,  uy  .  There  are  two  Interpretations  we  can  make 
for  the  u's.  In  terms  of  C^I,  we  argue  that  u's  monotonlcally  Increases  to¬ 
ward  one  as  the  C^I  Is  Improved,  I.e. ,  better  C^I  Improves  targeting  on  the 
average  and  hence  the  fire  tends  more  toward  aimed  fire  and  less  toward  area 
fire.  Another  Interpretation  is  that  of  an  exposed  force  attacking  a  concealed 
defender.  Suppose  X.  is  attacking  and  X  Is  defending.  Since  X  exposes  his  po¬ 
sition  during  attack,  X  may  a*™  Ms  fire  and  uy  tends  to  one.  On  the  other 
hand,  X  must  fire  Into  X's  defended  area,  not  knowing  exactly  where  X  Is,  so 
ux  tends  toward  zero.  This  Is  similar  to  a  case  studied  by  Brackney,  [1959]. 

In  an  Information  war,  the  state  variables  X  and  Y  represent  the  average 
Information  each  side  possesses  about  the  combat  environment.  We  are  inter¬ 
preting  Information  here  In  Its  theoretical  sense  as  lack  of  uncertainty. 
There  are  many  microscopic  informational  entitles  that  make  up  the  state  of  the 
combat  environment.  These  Include  your  own  and  the  enemies  force  locations, 
capabilities  and  even  his  Intentions.  The  true  state  of  the  environment  Is 
never  known  exactly  and  moreover,  without  continuing  intelligence  and  combat 
reports  to  update  or  replenish  the  Information,  it  will  decay  over  time  due  to 
the  dynamic  nature  of  the  environment. 

Each  Information  war  model  will  contain 

1)  A  natural  rate  of  loss  due  to  the  dynamic  properties  of  the  environ¬ 
ment. 

2)  Replenishment  of  Information  from  Intelligence,  sensor,  combat  and 
other  reports. 
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A  rate  of  loss  due  to  deliberate  acts  of  the  opposition  to  deceive, 
confuse,  jam,  etc.  This  Is  the  counter-03  terra  and  It  Is  In  this 
terra  that  the  three  Information  war  models  differ  In  their  assump¬ 
tions. 


X  ■  *b  X  -  p  (X  -  X)Y  +  Vv 

A  W  ^  A 

Y  -  -byY  -  px(Yq  -  Y)X  +  Vy 


(3a) 


X  -  -bxX  -  oyY  +  Vx 

•y  ,  .byY  -  axx  +  Vy 


(3b) 


X  *  -bxX  -  UyXY  +  Vx 

•  (3c) 

Y  «  -byY  -  pxXY  +  Vy 

In  all  three  models  bx  and  by  represent  natural  loss  rates  and  vx 
and  Vy  represent  replenishment  rates  for  X_  and  X  respectively. 

The  assumption  In  (3a)  Is  that  counter-c3  effectiveness  of  X  against  X 
depends  on  the  product  of  X's  knowledge  (Y)  and  X's  Ignorance  (X0  -  X). 
(X0  Is  the  maximum  possible  level  of  Information  of  X.)  The  premise  Is  that 
counter-C^  is  more  effective  against  an  already  confused  or  Ignorant  enemy 
than  against  an  Intelligent  one.  Furthermore,  that  one  Is  more  effective  In 
devising  counter-C^  the  better  Informed  he  Is  about  the  environment.  py 
and  px  are  parameters. 

In  (3b),  the  assumption  Is  that  the  enemies'  lack  of  knowledge  Is  Imma¬ 
terial  In  matters  of  counter-C^,  its  effectiveness  only  depends  on  one's  own 
knowledge  of  the  environment.  oy  and  ax  are  parameters,  here  assumed  to  be 
positive.  We  note  In  passing  that  for  negative  values  of  ax  and  ay,  (3b) 
are  Rlcharson's,  T1960]  liner  equations  for  arms  races. 
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In  (3c)  the  assumption  is  that  an  ignorant  or  already  confused  enemy  can¬ 
not  be  further  confused  by  counter-C^  acts,  but  that  in  fact,  the  effective¬ 
ness  of  such  tactics  actually  increases  with  his  as  well  as  one's  own  accurate 
state  of  knowledge  about  the  environment.  uy  and  yx  are  parameters,  pre¬ 
sumed  to  be  positive  but  possibly  greater  than  one.  Eq.'s  (3c)  were  those 
analyzed  for  equilibria  and  stability  in  the  paper  by  Moose,  [1980]  and  found 
to  be  environmentally  stable. 

The  final  model  we  wish  to  consider  is  in  the  area  of  intelligence  infor¬ 
mation  acquisition.  Like  the  other  information  models  X  and  Y  represent  the 
average  values  of  correct  information  about  an  uncooperative  environment.  The 
properties  to  be  modeled  are: 

1)  A  natural  rate  of  information  loss  due  to  the  dynamic  environment. 

2)  Replenishment  of  Information  through  normal  open  intelligence  gath¬ 
ering  channels. 

3)  Information  growth  through  uncooperative  channels  that  "leak"  more 
information  about  the  enemy  in  proportion  to  his  Ignorance  and  are 
better  "exploited"  by  one's  own  assets  in  proportion  to  one's  own 
knowledge.  The  latter  is  to  account  for  more  effective  targeting  of 
limited  acquisition  resources. 

X  *  -bxX  +  ux(Y0  -  Y)X  +  Vx 


(4) 


Y  «  -byY  +  Wy(X0  -  X)Y  ♦  Vy 

Eq.'s  (4)  are  the  intelligence  Information  evolution  equations  that  ac¬ 
count  for  these  proposed  properties.  ux  and  uy  are  parameters.  They  are 
presumed  to  be  positive  but  we  could  imagine  that  they  are  negative  If  the 
"leaks"  are  in  fact  "misinformation"  being  used  for  the  purpose  of  deception. 

In  summary,  the  five  systems  of  equations  represent  models  ofi 
Eq.  (2)  :  Mixed  attrition  combat. 
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Eq.  (3a)  :  Information  war;  counter-C^  depends  on  one's  own  know¬ 

ledge  and  opposition's  Ignorance, 

Eq.  (3b)  :  Information  war;  counter-C^  depends  on  one's  own  know¬ 

ledge, 

Eq.  (3c)  :  Information  war;  counter-C^  depends  on  the  product  of 
one's  own  knowledge  and  opposition's  knowledge, 

Eq.  (4)  :  Intelligence  information;  intelligence  grows  in  proportion 

to  enemy  ignorance  and  one's  own  knowledge. 
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Cannonical  Forms 


Eq.'s  (2),  (3a)  and  (3b)  are  of  a  cannonical  form  we  shall  call  "Type  1" 
equations.  Type  1  equations  are  designated  in  accordance  with  the  notation  of 
Eq.  (2),  as 

X  =  -bxX  ■  ay  Vy  Y  ■  “yf1  ■  U„)XY  +  Vx 

•  '  (5) 

Y  =  -byY  -  ax  Ux  x  -  «x(l  -  ux)XY  +  Vy 

The  parameters  bx  and  by  are  positive  and  have  the  interpretations  of  nat¬ 
ural  loss  rates.  Eqs.  (2)  are  obtained  from  (5)  by  restricting  ox  and  ay 
to  be  positive  and  ux  and  uy  to  be  greater  than  or  equal  to  zero  and  less 
than  or  equal  to  one.  Eq.  (3b)  Is  obtained  from  (5)  by  setting  yx  and  uy 
equal  to  one.  Eqs.  (5)  may  be  obtained  from  (3a)  by  the  transformations 


Since,  as  we  shall  explain  shortly,  X0  and  Y0  are  restricted  to  being 
greater  than  one,  uy  and  ,jx  will  be  greater  than  one  and  oy  and  <*x  will 
be  positive  in  this  transformation. 

In  summary,  Eq.  (5)  in  a  cannonical  form  for  Eq.  (2),  the  mixed  attrition 
model,  Eq.  (3a)  the  counter-C3  model  depending  on  the  enemies1  ignorance  for 
deceptive  successes,  and  Eq.  (3b),  the  counter-C3  model  that  is  Independent 
of  the  enemies'  state  of  knowledge.  The  range  of  the  parameter  pair  ux,  uy 
determines  the  model  to  which  (5)  applies.  Between  zero  and  one,  we  have  the 
mixed  attrition  model.  For  values  greater  than  or  equal  to  one  we  have  two  of 
the  counter-C3  models. 
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Eqs.  (3c)  are  of  a  unique  form  which  we  shall  designate  "Type  2"  equa¬ 
tions.  Again,  we  shall  Investigate  their  properties  in  the  two-dimensional 
parameter  space  yx,  uy  .  Eqs.  (4),  the  Intelligence  equations  are  also 
different  and  we  shall  designate  them  as  "Type  3"  equations. 


Stationary  Points 


In  general',  two-species,  quadratic  evolution  equations  may  have  four  sta- 

•  • 

tlonary  points,  that  is  points  in  the  X-Y  plane  for  which  X  and  Y  are  simulta¬ 
neously  zero.  However,  because  there  are  no  X*  or  Y?  terms  these  equations 
all  have  two  stationary,or  equll Ibriun,  points.  In  all  of  these  models,  the 
state  variables,  X  and  Y  only  have  Intelligent  Interpretation  when  they  are 
non-negative.  Furthermore,  the  purpose  of  the  replenishment  terms  is  to  make 
up  for  attrition  in  order  to  maintain  their  values  at  a  desired  level.  We 
shall  thus  presume  Vx  and  Vy  are  picked  to  establish  one  of  the  equilibrium 
points  In  the  positive  quadrant  of  the  state  space.  It  Is  mathematically  con¬ 
venient  to  pick  this  point  as  j  Xej,  Yej  |  by  suitable  scale 

changes  if  necessary.  Hereafter  we  shall  refer  to  this  as  the  unity  equilibri¬ 
um  point.  This  prescribes  a  fixed  relationship  between  the  model  parameters 
and  the  replenishment  terms.  Namely; 

Type  1: 


V  *  b  +  o 
x  x  y 


V  *  b  +  a 

y  y  x 


(7) 
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Type  2: 


V  *  b  +  u 
¥x  x  yy 

V  *  b  +  u 

y  y  wx 


Type  3: 

V  »  b  -  v„(Y  -  1) 

x  x  x  0  (g) 

Vy  3  by  -  uy(XQ  -  1). 

It  Is  easily  seen  that  (7),  (8)  and  (9)  substituted  Into  (5).  (3c)  and  (4) 
along  with  X  *  Y  =  1  make  X  and  Y  simultaneously  zero  for  the  Type  1.  Type  2 
and  Type  3  equations  respectively.  Given  that  unity  Is  an  equilibrium  point, 
one  may  find  the  other  equilibria  point  by  direct  algebraic  manipulation.  Let 

It  be  designated  |xe2,  Ye2^.  Then 
Type  1: 

(bu/av)(b>v  +  1)  -  Uv(Vax  +  1)  ^ 

x,2  •  Y  .*u  -  $  -  n.x^)d 


-  "x(bx/a»  *  11 

— »Ya  -  -  9 


Type  2: 


.  JL  (l  +  w  /b  ) 
ux  v  V  x' 

-  “  (1  +  vJK) 


give  the  locations  of  the  second  stationary  point  for  each  of  the  three  can- 
nonlcal  forms  In  terms  of  the  model  parameters. 


We  can  make  some  general  observations  about  the  locus  of  these  points. 
For  Type  1  systems,  the  second  equilibrium  point  may  lie  in  any  quadrant  of  the 
state  space.  However,  If  either  yx  or  yy  Is  greater  than  one.  It  cannot  be 
In  the  third  quadrant,  l.e.,  Xe2  and  Ye2  cannot  both  be  negative.  And 
for  px  >  1  +  by/ax  and  uy  >  1  +  bx/cty.  It  must  lie  In  the  first 
quadrant,  l.e.,  both  Xe2  and  will  be  positive. 

For  (3a),  the  first  counter-C3  model,  this  requirement  corresponds  to 
px  >  by  and  py  >  bx  according  to  the  transformation  required  by  (6). 
We  may  in  fact  use  (6)  to  define  the  second  equilibrium  point  of  Type  1  equa¬ 
tions  directly  in  terms  of  { X0,  Y0,  Px,  Py|,  the  counter-C3  parame¬ 
ters  of  (3a),  by  the  equivalent  form. 


+ 


(10a) 


For  large  values  of  px  and  py,  jxe2,  Yfi2j-»  —  »  Yo  "T —  )  *  For 

small  px  and  py,  the  equilibrium  point  moves  toward  plus  Infinity  until  at 
zero  the  equations  are  linear  for  which  there  is  just  one  equilibrium  point, 
unity.  In  between  ranges  of  px  and  py  can  move  the  equilibrium  point  into 
any  quadrant  with  the  behavior  quite  remarkable  near  the  singular  points. 

For  example,  as  py  -»  bx,  Xe2  -*  —  but  Ye2  -*  Y0  and  vice  ver¬ 
sa.  However,  If  py  -*  bx  and  px  -*  by  slmultaneoulsy,  then  Xe2  -*  X0 

and  Ye2  -*  Y0.  It  Is  for  values  of  px  and  py  near  these  singular 

points  that  we  find  j  Xe2,  Ye2|  outside  of  the  first  quadrant. 

In  (10),  the  original  form,  for  ux  and  uy  much  less  than  one,  a  likely 
range  for  (2)  the  mixed  attrition  Lanchester  equations,  the  second  equilibrium 
point  will  lie  In  the  fourth  quadrant,  and  be  Infeasible.  For  yx  -  py  *  1, 
the  case  of  (3b),  the  equations  are  linear  and  there  Is  only  one  equilibrium 
point,  unity.  The  locus  of  jxe2,  Ye2 1  Is  shown  In  Figure  1  for  several 
values  of  the  model  parameters,  that  might  correspond  to  (2)  and  In  Figure  2 
for  ranges  that  might  correspond  to  (3a).  It  Is  worth  mentioning  here  that 

X0  and  Y0  are  the  maximum  possible  levels  of  Information  X  and  Y  can  pos¬ 

sibly  possess.  If  unity  Is  to  be  an  equilibrium  point  that  Is  feasible,  then 
X0  and  Y0  must  be  greater  than  one,  a  condition  we  have  previously 

assumed. 

For  Type  2  systems,  Xe2  and  Ye£  are  always  negative.  Thus,  only 
one  equilibrium  point  Is  In  the  accessible  region  of  the  state  space.  The  neg¬ 
ative  point  Is  not  feasible  since  X  and  Y  cannot  be  negative  quantities.  Type 
2  equations  correspond  to  the  counter-C^  model  of  (3c)  vrfHch  postulates 
counter-C^  dependent  on  the  product  of  X/s  and  Y/s  knowledge.  This  was  the 
model  analyzed  by  Moose,  [1980]  where  the  same  condition  was  described. 
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The  second  equilibrium  point  for  Type  3  system  lies  In  the  third  quadrant 
providing  Vx  and  Vy,  the  constant  rates  for  replacing  lost  Intelligence  In¬ 
formation,  are  positive.  Even  If  Vx  and  Vy  are  zero;  Xe2.  Ye2f  * 

Jo.  »{•  1  1 

To  summarize:  When  their  first  equilibrium  point  Is  placed  at  unity.  In 
the  first  quadrant  of  the  state  space,  the  second  equilibrium  point  for  Type  2 
and  Type  3  systems  Is  strictly  negative,  l.e..  In  the  third  quadrant  of  the 
state  space.  The  second  equll  Ibrlum  point  of  Type  1  systems  may  be  In  any 
quadrant.  However,  for  sufficiently  large  values  of  ux  and  uy,  or  px  and 
py.  It  must  be  In  the  first  quadrant.  Only  the  first  quadrant  of  the  state 
space  is  considered  a  feasible  region  for  the  state  variables. 


Stability 


The  parameters  of  evolution  equations  determine  their  stability  near 
equll ibrlmn  points  as  well  as  the  location  of  the  equilibrium  points.  Stabil¬ 
ity  near  an  equilibrium  point,  called  neighborhood  stability,  refers  to  the 
dynamic  property  of  the  system  to  return  to  or  diverge  from  equilibrium  If  per¬ 
turbed  from  the  equll Ibrlun  position. 


In  order  to  investigate  this,  we  replace  the  state  variables  according  to 


X  »  8  +  X 
x  e 


Y  -  ay  +  Ye 


and  expand  (1)  about  |Xe,  Yej  In  a  Taylor  series.  Retaining  only  the  first 
order  terms  In  the  expansion,  which  Is  valid  for  small  8X  and  8y»  one  ob¬ 
tains  the  linear  state  equation 
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where  Ce,  called  the  conflict  matrix.  Is 


and  each  of  the  partial  derivations  are  to  be  evaluated  at  the  equilibrium 
point  under  Investigation.  [See  Moose,  [1983]  or  May,  [1974]  for  a  more  com¬ 
plete  description  of  the  perturbation  method  for  this  class  of  non-linear  evo¬ 
lution  equations.] 

The  conflict  matrix  determines  the  stability  or  Instability  of  the  system. 
If  all  its  eigenvalues  are  In  the  left  half  plane,  the  system  is  "stable"; 
I.e.,  perturbations  will  die  out.  If  one  or  both  of  Its  eigenvalues  are  In  the 
right  half  plane,  a  minor  perturbation  will  grow  and  the  system  will  diverge 
from  that  equilibrium.  In  this  case  It  is  said  to  be  "unstable."  If  Its  ei¬ 
genvalues  are  complex,  (they  must  be  complex  conjugates  If  they  are),  the  sys¬ 
tem  dynamics  will  be  oscillating.  The  oscillations  will  die  out  If  it  is  sta¬ 
ble,  but  will  grow  If  It  is  unstable.  If  the  eigenvalues  are  in  the  left  half 
plane  for  all  possible  values  of  the  parameters,  the  system  Is  said  to  be  "en¬ 
vironmentally  stable."  If  for  some  values  of  the  parameters,  the  system  Is 
stable  but  for  others  It  Is  unstable,  the  system  is  said  to  "environmentally 
unstabl  e. " 

The  conflict  matrices  for  the  three  cannonical  forms  are; 

Type  1: 

"  bX  ■  “y^1  "  My)Ye  “ay  »y  "  «y(1  '  Wy)Xe 

Ce  "  d«> 

*  <*x  ux  -  ax(l  -  ux)Yfi  -b  -  ax(l  -  wx)Xe 
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Type  1  (alternate  form); 


-b  +  p  Y 
x  wy  e 


VYo  "  V 


VXo  -  Xe> 


"by  +  px  Xe 


Type  2 


-b  -  u  Y 
x  My  e 


"Mx  Ye 


'wy  Xe 


-by  -  wx  Xe 


Type  3 


-bx  *  “x(¥o  -  V 


“e 


-by  *  ,y(X#  -  Xe) 


Equations  (16)  through  (18)  show  how  the  conflict  matrix  and  hence  system 
stability  depends  on  the  parameters  of  each  of  the  cannonlcal  types.  Recall 
that  jxei,  Yej  }  *  |l,  l|  Is  an  equilibrium  point  located  In  the  feasi¬ 
ble  region  of  the  state  space.  In  order  to  obtain  the  conflict  matrices  near 
this  point,  we  simply  put  Xe  *  Ye  =  1  In  (16)  through  (18),  and  Investigate 
the  properties  of  the  eigenvalues.  However,  j  Xe2,  Ye2  j  Is  also  an 
equilibrium  point  and  we  must  substitute  the  appropriate  equations  from  the 
previous  section  Into  (16)  through  (18)  In  order  to  determine  stability  there. 

If  equilibrium  point  two  as  well  as  the  unity  equilibrium  Is  In  the  feasi¬ 
ble  region,  as  Is  frequently  the  case  In  Type  1  systems,  then  there  Is  the  pos¬ 
sibility  of  two  feasible  stable  steady  states.  Or  one  or  the  other  may  be  sta¬ 
ble  and  the  other  unstable,  In  which  case  the  system  will  seek  the  stable  solu¬ 
tion.  For  Type  2  and  3  systems,  the  stability  of  the  unity  point  Is  most 
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Important  because  the  other  point  Is  not  feasible.  However,  if  the  unity  point 
Is  unstable,  and  the  unfeasible  point  Is  stable,  then  we  can  forecast  that  both 
the  state  variables  will  die  out,  that  Is,  tend  toward  negative  values.  Instead 
of  growing  without  bound. 

The  eigenvalues  of£e  are  the  roots  of  the  characteristic  polynomial 

D(p)  «  p2  -  (Cn  +  C22)p  +  (Cn  C22  -  C12  C21)  (19) 

where  the  coefficients  are  the  elements  of  Ce. 

The  roots  of  (19) 


are  functions  of  the  system  parameters  according  to  (16)  through  (18).  We  note 


that  If 

(Cn  +  c22)  <  0 

and 

C11  C22  >  C12  C21 


(21) 


then  Re[pj  ]  and  Re[p2]  are  less  than  zero,  l.e.,  pj  and  pj  are  In 
the  left  half  plane,  and  the  system  Is  stable.  Eq.  (21),  which  is  called  the 
Routh  Test  [see,  e.g.,  Truxall,  [1972]]  Is  a  necessary  and  sufficient  condition 
for  stability.  We  also  note  from  (20)  that  Cjj  C^i  >  0  is  a  sufficient 
condition  for  both  of  the  roots  to  be  real. 


Stability  of  Type  1  Systems  at  Unity  Equilibrium 

Type  1  systems  are  stable  at  unity  equilibrium  for  all  ux  and  py  less 
than  one  when 

twx  -  (by/ox  +  1)]  [jiy  -  (bx/ay  +1)]>1.  (22) 
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Eq.  (22)  Is  the  second  condition  of  (21),  since  the  first  Is  always  satisfied. 
The  stable  and  unstable  regions  determined  by  (22)  are  Illustrated  In  Figure  3. 
Note  that  when  natural  loss  rates  are  equal  to  or  greater  than  the  external  at¬ 
trition  coefficients,  b x/<*y  >  1  and  by/ax  >  1,  the  system  Is  stable  for 

all  positive  values  of  jix  and  uy  less  than  one.  On  the  other  hand.  If 
bx/ay  and  by/ax  become  very  small,  then  the  unity  equilibrium  point  Is 

unstable  for  all  values  of  ux  greater  than  bx/ay  and  uy  greater  than 

by/ax.  In  Figures  3a,  3b  and  3c,  we  have  sketched  the  left  side  of  (22) 
for  three  specific  values  of  bx/ay  and  by/ax  In  order  to  Illustrate 

this  point  graphically. 

The  quantities  b/a  are  related  to  the  efficiency  of  warfare.  When  warfare 
Is  very  efficient  b/a  Is  small  for  both  sides  and  war  Is  fundamentally  unstable 
at  Its  unity  equilibrium  point;  that  Is,  the  force  levels  cannot  stay  at  their 
planned  military  levels.  What  happens  with  larger  values  for  v,  corresponding 
to  greater  efficiency  Is  that  It  Is  unstable  for  relatively  larger  values 
of  b/a  corresponding  to  lower  levels  of  warfare  efficiency. 

Since  the  combat  force  levels  cannot  remain  at  their  planned  levels,  when 
the  system  Is  unstable,  we  must  ask,  what  will  the  system  trajectory  be?  This 
brings  the  stability  properties  of  the  second  equilibrium  point  Into  play. 
Because  Xe2,  Ye2,  which  depend  on  the  system  parameters,  appear  explic¬ 
itly  In  the  conflict  matrix.  Its  eigenvalues  are  much  more  complicated  func¬ 

tions  of  the  parameters  than  the  eigenvalues  at  unity  equilibrium.  A  discus¬ 
sion  of  these  properties  will  be  deferred  to  a  separate  paper. 

For  ux  and  uy  greater  than  one,  Type  1  equations  represent  the 

counter-C^  model  of  Information  war  of  (3a)  and  we  can  determine  the 


IQ 


conditions  from  (16a)  and  (21)  for  stability  to  be 


o  +  0  (.  b  +  b 
px  py  x  y 

and  v  /  v  (23) 

fe  '  ‘)  fe  ‘ ')  >  (x°  ■  1)(Y°  ‘ 11 

It  Is  easy  to  show  that  the  second  condition  Implies  the  first  so  that  the 

second  condition  Is  both  necessary  and  sufficient  for  the  system  to  be  stable. 
As  shown  In  Figures  4a  and  4b,  this  occurs  for  values  of  px  less  than  by 
and  py  less  than  bx.  As  with  the  Lanchester  model  of  mixed  attrition,  this 
Information  war  model  becomes  unstable  at  unity  equilibrium  trfien  the  counter¬ 
ed  effectiveness  jpx,  py|  exceeds  the  natural  loss  rates,  jbx,  by  . 

Examples  of  stable  and  unstable  phase  trajectories  are  shown  In  Figures 

5a,  5b,  5c  and  5d  for  Type  1  systems.  Type  1  systems  are  environmentally 
unstable  systems  as  defined  earl ler. 

Finally,  we  note  that 

c12  C21  '  V  »x  0y^xo  '  We  '  »  >  0  (24) 

for  Type  1  systems  at  unity  equilibrium.  Therefore  the  eigenvalues  are  real 
numbers  and  the  trajectories  are  non-oscillating  near  there. 


so  that  the  system  Is  stable  when 


S*  •  /  <  1  <26> 

y  x 

l.e.,  when  the  warfare  efficiencies  are  relatively  low.  When  efficiencies 
become  too  great,  the  system  becomes  unstable  and  diverges  toward  vrfiichever 
side  gains  the  Initial  advantage.  There  is  no  second  equilibrium  point.  The 
regions  defined  by  (26)  are  shown  In  Figure  6  for  the  jax,  oy  |  parameter 
space. 


Stability  of  Type  2  Systems 

When  we  apply  the  Routh  Test,  (21),  to  the  conflict  matrix  for  Type  2  sys¬ 
tems,  (17),  at  unity  equilibrium,  we  find  that  the  system  is  environmentally 
stable  at  unity  equilibrium  since 


and 


(bx  +  by  +  ux  +  uy}  <  0 


bb+ub+ub>0 
x  y  px  x  My  y 


(27) 


for  all  possible  values  of  the  parameters. 

At  equilibrium  point  two,  which  is  always  in  the  third  quadrant,  we  find, 
with  the  aid  of  (11)  that  it  is  always  unstable  since 


C11  +  C22 


wx  bx^by  +  uy  *Vbx  >  0 


and  the  first  condition  of  (21)  is  always  violated.  Furthermore 


C12  C21 


ux  Uy  Xe  Ye  >  0 


(28) 


(29) 


so  the  solutions  of  the  linear  state  equations  are  non-osclllatory. 

To  summarize;  Type  2  systems,  the  counter-C3  model  (3c)  and  the  model 
studied  by  Moose,  [1980]  for  counter-C3,  are  environmentally  stable  at  unity 
equilibrium,  and  non-osclllatory.  ihelr  other  equilibrium  point  which  Is  In 
the  third  quadrant  and  is  therefore  unfeasible,  is  always  unstable.  Thus,  we 
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may  anticipate  Type  2  systems  to  always  tend  back  toward  unity  equilibrium, 
even  from  quite  far  away.  In  fact,  they  are  probably  globally  stable  at 
unity. 


Stability  of  Type  3  Systems  at  Unity  Equlllbri 


urn 


At  unity  equilibrium.  Type  3  systems  are  stable  if 


“x<Yo  -  »  ♦  »y(X0  -  ')  <  bx  *  b. 


and 


b 

“  -  (Y  -  1) 


-  (Xft 

u  o 


i) 


<  l 


(30) 


Eq.  (30)  is  quite  similar  to  (23),  the  stability  conditions  at  unity  for 
Type  la  systems,  but  the  second  inequality  is  reversed.  Figure  7  illustrates 
how  the  conditions  of  (30)  generate  a  region  of  stability  which  is  a  diagonal 
strip  in  the  |ux,  uy  |  parameter  space.  For  either  too  small  or  too  large 
values  of  ux  and  uy,  the  system  is  unstable. 

We  know  that  the  second  equilibrium  point  of  Type  3  systems  always  lies  in 
the  third  quadrant  and  is  therefore  not  feasible.  It,  like  unity  equilibrium, 
is  environmentally  unstable.  Its  properties  will  also  be  described  more  fully 
in  a  subsequent  report.  However,  we  note  here  that 


C12  C21 


u  X  Y  >  0 
x  y  e  e 


(31) 


holds  for  both  points  so  that  Type  3,  like  Type  2  systems,  are  non-oscillatory. 

A  summary  of  the  principal  system  properties  for  each  of  the  five  models 
is  presented  below  in  Table  1. 
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Type  1 

Type  2 

Type  3 

3 

Mixed  Attrition  Counter-C 
E<^(2) _ Ed.  (3a) 

3 

Counter-C 
Eq.  (3b) 

3 

Counter-C 
Eq.  (3c) 

Intell igence 
Eq.  (4) 

1) 

Unity  Equil ibrium 
Stabil  ity 

E.U. 

E.U. 

E.U. 

S. 

E.U. 

2) 

Equil  ibrium  Point 
Two  Locus 

I, II, III, IV 

I, II, III, IV 

n/a 

Ill 

Ill 

3) 

Equilibrium  Point 
Two  Stability 

E.U. 

E.U. 

n/a 

u. 

E.U. 

4) 

Oscill atory 

Sol utions 

yes 

(Equl ibm  2) 

yes 

(Equl ibm  2) 

no 

no 

no 

Table  1 


Summary  of  Principal  System  Properties 


(E.U.  =  Environmentally  Unstable,  U  *  Always 
Unstable,  S  =  Always  Stable,  State  Space  Quad¬ 
rants  I,  II,  III,  IV.) 


Discussion 


The  five  models  we  have  Investigated  in  this  paper  are  all  generalized 
Lotka-Vol terra  equations  for  two-species  systems.  In  spite  of  these  restric¬ 
tions,  their  interpretation  in  the  C3I  context  seems  meaningful  and  the  na¬ 
ture  of  their  dynamic  behavior  quite  remarkable. 

The  mixed  attrition  Lanchester  type  model  with  constant  replenishment  is 
quite  useful  to  evaluate  the  Increased  effectiveness  of  forces  as  they  become 
more  precisely  targeted  against  the  enemies'  forces.  The  Interpretaton  of  the 
parameter  py  for  side  Y  can  be:  "uy  increases  from  zero  toward  one  as  Y's 
targeting  accuracy  improves  with  better  C3I,"  or  "py  decreases  toward  zero 
as  counter-C3  techniques  by  X  confuse  or  decoy  Y's  aimed  fire"  or  "py  is 
large  if  X,  in  order  to  attack,  must  expose  his  positions  to  Y  who  has  accurate 
targetable  weapons."  (For  example,  mines  are  not  targetable  weapons  but  laser- 
aimed  anti-tank  missiles  are.)  All  three  of  these  Interpretations  fcr  the  p 
parameter  in  the  mixed  attrition  model  are  C3I  or  counter-C3  related.  We 
also  note  that  increasing  p  tends  to  make  the  system  unstable.  One  would  sus¬ 
pect  that  the  side  with  the  greatest  p  should  win  vrfien  instability  occurs.  Al¬ 
though  this  may  happen.  It  is  not  guaranteed  and  in  fact  what  does  transpire 
depends  critically  on  the  location  and  stability  of  the  second  equilibrium 
point  and  on  the  initial  conditions  at  the  onset  of  the  conflict.  We  will  In¬ 
vestigate  this  phenomenon  In  more  detail  in  a  separate  report. 

The  three  different  models  for  information  war  Involving  counter-C3 
measures  show  the  critical  Importance  of  our  assumptions  about  how  the  counter- 
C3  techniques  actually  effect  the  opposition's  efforts  to  keep  his  Informa¬ 
tion  current  and  accurate.  If,  as  was  assumed  In  the  previous  paper  [Moose, 
[1980]],  these  activities  are  effective  against  a  well  Informed  enemy,  but 


ineffective  against  an  ignorant,  or  already  confused  one,  then  the  dynamics  are 
those  of  the  Type  2  systems.  These  dynamics  are  remarkably  simple  and  essen¬ 
tially  Independent  of  the  actual  values  of  the  model  parameters  insofar  as  the 
gross  system  behavior  is  concerned.  On  the  other  hand,  if  counter-C3  methods 
are  in  fact  more  effective  against  an  enemy  as  he  becomes  more  confused  or  ig¬ 
norant,  such  as  might  occur  through  input  overloading  of  his  sensors  and  proc¬ 
essing  equipments,  then  the  dynamics  are  those  of  Type  la  systems.  These  dy¬ 
namics  have  two  major  modes,  stable  and  unstable,  and  the  mode  of  behavior  is 
critically  dependent  on  the  actual  model  parameters. 

Finally,  in  our  intelligence  Information  model,  the  behavior  corresponds 
to  Type  3  systems.  Here  again  there  are  two  modes  of  behavior  but  only  one 
possible  stable  point  in  the  feasible  region.  Therefore,  either  both  sides  are 
maintained  with  a  fairly  constant  level  of  knowledge  about  the  other  or  one 
side  tends  to  know  a  great  deal  while  the  other  knows  very  little.  The  first 
mode  prevails  when  intelligence  comes  primarily  from  analysis  of  non- restricted 
or  uncontrollable  sources  of  Information  about  the  enemy.  The  second  condition 
prevails  when  Intelligence  comes  more  predominately  from  breaks  in  the  enemy's 
security;  through  espionage,  communication  intercepts,  etc.  This  model  too 
needs  further  examination  as  the  stability  of  its  second  equilibrium  point, 
which  is  in  the  third  (non-feaslble)  quadrant,  still  has  a  great  Influence  of 
the  actual  trajectory  of  this  type  of  information  war. 
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figure  1 

Locus  of  Equilibrium  Pt.  Two 
Mixed  Attrition  Type  1  Systems 


11 


FIGURE  2 

Locus  of  Equilibrium  Pt.  Two 
Information  War  Type  1  Systems 


FIGURE  3a 

Stability  Regions  for  Unity  Equilibrium 
Mixed  Attrition  Type  1  Systems 


FIGURE  3b 

Stability  Regions  for  Unity  Equilibrium 
Mixed  Attrition  Type  1  Systems 


FIGURE  3c 

Stability  Regons  for  Unity  Equilibrium 
Mixed  Attrition  Type  1  Systems 


FIGURE  4a 

Stability  Regions  for  Unity  Equilibrium 
Information  War  Type  1  Systems 


FIGURE  4b 

Stability  Regions  for  Unity  Equilibrium 
Information  War  Type  1  Systems 
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FIGURE  5a 

Stable  Phase  Trajectory 
Mixed  Attrition  Type  1  Systems 


FIGURE  5b 

Unstable  Phase  Trajectory 
Mixed  Attrition  Type  1  Systems 


FIGURE  5c 

Stable  Phase  Trajectory 
Information  War  Type  1  Systems 


FIGURE  7 

Stability  Regions  for  Unity  Equilibrium 
Information  War  Type  3  Systems 


DISTRIBUTION 


No. 


1.  Library,  Code  0212 
Naval  Postgraduate  School 
Monterey,  California  93943 

2.  Dean  of  Research,  Code  012 
Naval  Postgraduate  School 
Monterey,  California  93943 

3.  Director  of  NET  Assessment 
Attn:  LtCol  Fred  Giessler 
Office  of  Secretary  of  Defense 
Room  3A930,  Pentagon 
Washington,  DC  20301 

4.  Dr.  Joel  Lawson 
Technical  Director 

Naval  Electronics  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20360 

5.  Office  of  Chief  of  Naval  Operations 
Attn:  Dr.  Daniel  Schutzer,  OP  009T 
The  Pentagon 

Washington,  OC  20301 

6.  Professor  K.E.  Woehler,  Code  61Wh 
Department  of  Physics  &  Chemistry 
Naval  Postgraduate  School 
Monterey,  California  93943 

7.  Professor  Paul  H.  Moose,  Code  62Me 
Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93943 

8.  Dr.  Thomas  P.  Rona 
Boeing  Aerospace  Co. 

P.0.  Box  3999 

Seattle,  Washington  98124 

9.  Office  of  Naval  Research 
Attn:  Dr.  Charles  Hoi  and 
Arlington,  Virginia  22217 

10.  Laboratory  for  Information  and  Decision  Systems 
MIT  Attn:  Dr.  Robert  R.  Tenney 

Cambridge,  Massachusetts  02139 


/ 

t 


Copies 

4 

1 

1 

1 

1 

1 

25 

1 

1 

1 


40 


